Theory of helicity-sensitive terahertz radiation detection by field effect transistors 
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Within the two antenna model, we develop a theory of the recently observed helicity-sensitive 
detection of terahertz radiation by FETs. The effect arises because of the mixing of the ac signals 
produced in the channel by the two antennas. We calculate the helicity-dependent part of the 
photoresponse and its dependence on the antenna impedance, gate length, and gate voltage. 



Introduction. Field effect transistors (FETs) can be 
used for efficient detection of terahertz radiation 
The standard model [l[ assumes that the radiation is 
coupled to the transistor by an effective antenna which 
generates an ac voltage predominantly on one side of the 
transistor (e.g. between source and gate contacts). Since 
the ac gate-to-channel voltage modulates both the elec- 
tric field and the electron concentration in the channel, 
the current density will contain a dc component, which 
leads to a photoresponse in the form of a dc source- 
drain voltage AU proportional to the radiation intensity. 
Terahertz radiation detection and imaging by FETs was 
demonstrated in many experimental works, see reviews 
in Refs. Hi]. 

For a number of applications it is essential to charac- 
terize the polarization of terahertz radiation. The strong 
sensitivity of FET detectors to linear polarization has 
been demonstrated [H-i] and was shown to be caused 
by the anisotropic sensitivity of the effective antennas 
(formed by bonding wires and contact pads). However 
the problem of general terahertz polarimetry, including 
the determination of circular or elliptical polarization is 
still open. The first step in this direction was made in 
the recent experiments of Drexler et al [9( , who discovered 
the strong sensitivity of the FET detector to the helicity 
of terahertz radiation, i.e. to the sign of the phase shift 
between the x and y components of the radiation field. 

The qualitative explanation of this finding is based on 
the assumption (which was verified experimentally) that 
there are two effective antennas connected to the source 
and the drain sides of the transistor. One of them is pre- 
dominantly sensitive to the x-polarization and the other 
one - to y-polarization. For circularly polarized radia- 
tion, the ac voltages of these antennas are thus phase- 
shifted by 90°. To detect the helicity, the photoresponse 
must depend on this phase shift. This in turn is possible 
only because of the mixing of the ac signals produced in 
the channel by the two antennas. 

To understand the conditions for such mixing one 
should take into account that (i) the condition lot <C 1 
is normally satisfied, where ui is the radiation frequency 
and r is the momentum relaxation time and (ii) under 
this condition the ac current injected at the source leaks 
to the gate through the distributed gate-to-channel ca- 
pacitance on a distance defined by the ac leakage length I 
[3| . Thus the mixing of ac signals produced by the source 
and drain antennas can occur only if the gate length L is 
comparable to the leakage length I. If L 3> I, the mixing 



becomes exponentially small. 

This article presents a theoretical study of the helicity- 
sensitive detection of terahertz radiation by FETs within 
the two-antenna model 0] explained above. We show 
that to describe the required mixing of the two signals one 
has to abandon the usual assumption that the antennas 
can be considered as ac voltage sources (i. e. having zero 
ac impedance). Such mixing can occur only when the 
antennas impedances are comparable to or greater than 
the transistor impedance, and in particular when the 
antennas are current sources (infinite impedance). We 
calculate the FET photoresponse, including its helicity- 
sensitive component, as a function of the gate length L 
and gate voltage V g . This component has a maximum 
for certain values of these parameters. 
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FIG. 1: Schematics of the high-frequency circuit of FET with 
source and drain antennas. 



The problem. We consider a FET with antennas at 
both the source and the drain sides, as presented in 
Fig. 1. The antennas are assumed to generate ac volt- 
ages U a cos oot and £/&cos(wi + 0) respectively, 9 being 
the phase shift. To simplify the calculations we assume 
that the antennas have equal impedances at the radia- 
tion frequency lo denoted by Z. We also assume that a 
dc voltage V g is applied between the source and the gate, 
while the dc condition at the drain is open circuit. Wc 
consider the case oot <C 1. 

The basic equations are the continuity equation and 
the Ohm's law 0: 



dp 
dt 



dx 



= 0, 



8U 

dx ' 



(1) 



where p and j are the electron charge and current densi- 
ties in the channel, a = pp is the 2D channel conductiv- 
ity, fi is the electron mobility, U = V g — Vth is the gate 
voltage swing, V t h is the threshold voltage. 

We assume that the spatial variation of U is large com- 
pared to the gate-to-channel separation. Then the values 
of p and a are determined only by the local value of U 
(the gradual channel approximation). For an open chan- 
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nel (U > 0), the plane capacitor formula is applicable: 

p = CU. (2) 

Here C is the gate-to-channel capacitance per unit area. 
This relation is not valid below threshold (U < 0) as well 
as in the vicinity of the point U = 0. To simplify the 
following calculations we will derive our results assuming 
Eq. @ to hold and the electron mobility p to be constant. 
We will then present the modifications needed to account 
for the general dependences p{U) and a(U). 

Combining Eqs. (Q]) and @ we obtain the non- linear 
diffusion equation 



dU 
~dt 



p d 2 U 2 
2 dx 2 ' 



(3) 



We expand U(x, t) = U^ + U\(x, t) + U 2 (x, t) up to the 
second order of perturbations induced by radiation, i.e. 
Ui is the hrst order correction to Uq and U 2 is the second 
order correction. In the first order Eq. @ yields: 



(4) 



The boundary conditions corresponding to Fig. 1 are: 



U a cos tut — j Z = Ui 
U b cos(ujt + 6) +jZ ■ 



at x = 0, 
at x == L. 



(5) 



where Z is the antenna impedance. 

From now on we will use dimensionless variables, as- 
suming that the voltages are in units of Uq, the coordi- 
nate x and th e gate l ength L are in units of the leakage 
length I = \/2a /luC = ^2p,Uo/ui [3(, and time t is in 
units of l/u. In these units Eqs. ((!]), ([5]) are rewritten 
as: 



dUi _ l d 2 U! 
dt ~ 2 dx 2 ' 
C/aCosi + Cfr = Ux atx = 0, 
U b cos(t + 9)- = Ui &tx = L, 



(G) 
(7) 



where £ = Zctq/I is the dimensionless antenna 
impedance, <to is the channel conductivity in the absence 
of radiation. In fact, £ is the ratio of Z and the resistance 
of the rectifying part of the channel of length I. 

To find the photoresponse voltage we need to consider 
the second order equation following from Eq. ([3]): 



dU 2 
dt 



d 2 U 2 
dx 2 



1 d 2 u 2 

2 dx 2 



(8) 



We average this equation over time and integrate twice 
over x [|| . The constant that appears after the first inte- 
gration is the dc current through FET which should be 
zero due to the dc open circuit condition on the drain 
side. Thus the photoresponse is given by 



AU = (U 2 (L) - U 2 (0)) 



1 r 
2 



{um 2 )-(Ui{l? 



(9) 



where the angular brackets denote time-averaging. If the 
impedance of antennas is negligible (i.e. Z = 0), then 
Ui(x — 0) = U a cos cot and U\(x = L) = Ubcos(ojt + 9). 
For this case we obtain the simple result: 



AU= 1 -(U 2 -U 2 ). 



(10) 



Therefore, if the antennas are ac voltage sources, their 
signals are not mixed, and consequently there is no sen- 
sitivity to the phase shift 9. 

Solution. The general solution of equation Eq. (|6|) is 



Ui(x,t) = — 



A e (l+i)(*-L/2) + Be (l+i)(L/2-x) 

where c.c. means complex conjugated. Therefore 



(Ui(x,t) 2 ) 



Ae (l+t)(x-L/2) + Be (l+i)(L/2-x) 



C.C, 
(11) 

■ (12) 



The boundary conditions Eq. ([7]) give the equations for 
the coefficients A and B: 



aA + /3B = U a , 
/3A + aB = U b e if> , 



(13) 



with a = (1 - rj)e- x and /3 = (1 + i])e x , A = (1 + i)L/2, 
i] = (I + Using Eq. (fT2j) we can now present the 
photoresponse voltage given by Eq. ((9]) in the form 



AU 



Ae 



Be> 



\Ae x + Be' 



Re 



(B - A) (A* + B*) shA chA* , (14) 



where the star denotes complex conjugation. With the 
help of Eq. (|I3p we obtain: 



AU = Re 



U a -U b e w fU a + U b e 



shA chA* 



(15) 



\3 -a \ a + P 
Finally, the expression for the photoresponse AU reads 



AU 



-^{Ul-UD + ^UaUbSmg, (16) 



where the coefficients Fq and F\ are determined by the 
dimensionless antenna impedance £ and the dimension- 
less gate length L, entering via the parameters rj and A: 



F = Re(G), 
G 



Fi = Im(G), 



(1 + ??cthA)(l +?/thA)* 



(17) 



We note that the dependence of AU on U a , U b , and 9 
in Eq. (|16|) naturally follows from the symmetry of our 
problem resulting from the simplifying assumption that 
the impedances of source and drain antennas are equal. 
In this case, interchange of U a and U b together with the 
reversal of the phase shift 9 should obviously change the 
sign of the photoresponse AU. 
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In particular, the interference term containing U a Ub, 
changes sign when 9 is replaced by —0. We recall that 
source and drain antennas are assumed to predominantly 
respond to x- and y-polarizations respectively, meaning 
that the contribution of the interference term has oppo- 
site signs for right and left circular polarization. If the 
impedances of the antennas were different, sin 9 would be 
replaced by sin(# + #o), where 9q is an additional intrinsic 
phase shift related to the impedance difference. 

Let us consider some limiting cases. 

a) Antennas are ac voltage sources, i.e. Z = and 
r] = 0. Then Eqs. $Tjo$) and (HTJ) give F = 1, Fi = and 
we retrieve the simple Eq. (|10[) . 

b) Antennas are ac current sources, i.e. £ 3> 1. In 
this case it is convenient to introduce the antennas' ac 
current amplitudes j a = U a /\i]\ and ji, = Ub/\r/\. Using 
Eqs. ([HI) and CLZ]), we obtain: 



A(7 = ~f (j% - f b ) + ~/ij ji,sinf 



where 



fo = 



sh 2 L — sin 2 L 
sh 2 L + sin 2 L' 



fi = 



2shL sin L 
sh 2 L + sin 2 L 



(18) 



(19) 



c) Long gate, L > 1. In this case the interference 
between source and drain ac signals vanishes, and we 
obtain the result similar to Eq. (fTO)) : 



AU 



4 |1 + 77| 2 ' 
which accounts for the finite antenna impedance. 




FIG. 2: The coefficient F\ defining the helicity-sensitive part 
of the photoresponse AU as a function of the gate length 
L (in units of leakage length I) for different values of the 
dimensionless antenna impedance £, as indicated. For conve- 
nience of presentation, F\ is multiplied by the factor 1 1 + »7| 2 = 
|1 + (1 + i)C\ 2 - Note the change of sign around L — ty. Inset: 
dependence of Fi on f at L = 1. 

We are mostly interested in the helicity-sensitive part 
of the photoresponse described by the coefficient F\ in the 



second term of Eq. (16). Fig. 2 presents F\ as a function 
of the gate length for different values of the dimensionless 
antenna impedance Q. Maximum sensitivity to radiation 
hclicity occurs when the gate length is around the leakage 
length I. As a function of ( for a given L, the coefficient 
F\ has a maximum around £ = 1, which corresponds to 
ac impedance matching, see the inset in Fig. 2. 

Modifications for sub-threshold regime. So far we con- 
sidered detection by an open channel (U — V g — Vth > 0) 
when Eq. 2 is valid. However experimentally the maxi- 
mum of the photoresponse occurs below threshold 0, [E[ . 
Fortunately, our theory can be easily extended to arbi- 
trary gate voltages including sub-threshold values. 

Our solution of the basic Eq. 1 relied on the simple 
expressions p = CU and a = pCU for the charge density 
and conductivity respectively, as well as on the implicit 
assumption that the mobility y, does not depend on U. 
This is a good approximation in the open channel regime 
(U > 0). Generally, one has to consider p and a as 
non-linear functions of U not neccessarily proportional 
to each other. Then the derivatives dp/dt and da/dx, 
appearing in Eq. 1, should be expessed as 



dp dp dU 
~dt = dU dt ' 



da 

dx 



da_dU_ 
dU~dx~' 



(20) 



For the simplest situation (one ac voltage source an- 
tenna exciting a long gate FET) , this consideration leads 
to an important generalization of the expression AU = 
U% /4C/ El , which reads 0, El : 



(21) 



This formula which replaces 1/JTo by [dhxer J dU)u—u is 
valid for arbitrary gate voltage, it also takes into account 
the possible dependence p,(U). 

Using Eqs. 20, it can be shown that in our problem the 
generalisation to arbitrary dependences p(U) and o~(U) 
can be reduced to the simple change in units of length 
(leakage length I) and voltage (gate voltage swing Uq)- 



/ w 1 dp \ - 1 / 2 
V 2p, p dU J u=u ' 



U 



--Y 1 ■ (22) 
adUJu=u y ' 



With these modifications of units, the previous results in 
Eqs. (16, 17) remain valid. 

Gate voltage dependence. Now we can analyze the de- 
pendence of the helicity-sensitive part of the photore- 
sponse on the gate voltage swing Uq. This dependence is 
governed by the interplay of several factors: (i)change of 
the ac impedance of the active parts of the transitor and 
corresponding variation of the parameter £, (ii) change of 
the ratio L/l due to the variation of the leakage length I, 
and (iii)variation of the voltage unit in Eq. 22. 

Assuming the mobility p to be constant, we will use 
for p(U) the semi-empirical formula (l2| : 



p = CU* ln[l + exp(U/U*)], 



(23) 



4 



where U* — akT/e, T is the temperature, and a is a 
phenomenological parameter on the order of 1. This for- 
mula reduces to Eq. 2 for positive U 3> U* and gives an 
exponential decrease of p for negative U. 

As explained above, the generalized expression for the 
photoresponse is given by Eq. 16 provided the voltages 
U a , Ub, and AU are measured in voltage units in Eq. 
22. However this unit is inconvenient because it itself 
depends on Uq. For this reason we will switch to a con- 
stant voltage unit U*. Then Eq. 16 will have the same 
form, except that the coefficients Fq and Fi, should be 
replaced by 



F U*(^ 



a dU I u=u ' 



Tx = FxU 



KadU 



J u=u ' 



We have performed numerical calculations of T\ using 
Eq. 23. Fig. 3 shows the dependence of J-i on the dimen- 
sionless gate voltage swing Uq/U* for several values of the 
antenna impedance. Since the dimensionless impedance 
£ changes with gate voltage, we label the curves by the 
values of £ at threshold (Uq = 0). Also, the gate length 
L was taken equal to the leakage length I at threshold. 

One can see that the hclicity-dcpcndcnt contribution 
to the photoresponse has a maximum around or below 
threshold. The decrease at the left side is due to the 
growth of the channel resistance. As a consequence, the 
dimensionless antenna impedance decreases approaching 
the limit Q = with no sensitivity to helicity, see the case 
a) above. The decrease at the right side results from the 
drop of the ratio L/l (see Fig. 2) caused by the increase 
of the leakage length I. 

Consistent with our results, an increase of the helicity- 
dependent contribution to the FET photoresponse when 
moving towards negative gate bias was observed experi- 



mentally However, the existing experimental data is 
not sufficient for a detailed comparison with our results. 
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FIG. 3: The helicity-sensitive contribution Ti as a function of 
Uo/U* for different values of as indicated, and L — I (both 
£ and I taken at Uo = 0). 

In summary, we have presented a theory of helicity- 
dependent terahertz detection by FET based on the 
model with two antennas connected to the source and 
drain sides of the transistor and sensitive to orthogonal 
linear polarizations. We have found the dependence of 
the helicity-sensitive component of the photoresponse on 
the antenna impedance, gate length, and gate voltage. 
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